Spatial and spectral shape of 
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We develop a mean-field theory of the spatial profile and the spectral properties of polariton 
condensates in nonresonantly pumped semiconductor microcavities in the strong coupling regime. 
Predictions are obtained for both the continuous-wave and the pulsed excitation regimes and the 
specific signatures of the non-equilibrium character of the condensation process are pointed out. 
A striking sensitivity of the condensate shape on the optical pump spot size is demonstrated by 
analytical and numerical calculations, in good quantitative agreement with recent experimental 
observations. 
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First evidences of Bose-Einstein condensation (BEC) 
in a solid-state system have been recently reported in a 
gas of exciton-polaritons in a semiconductor microcavity 
in the strong coupling regime P, 0, H, 0) S Hi ■ In addi- 
tion of being a remarkable example of an exciton conden- 
sate, this system opens interesting perspectives towards 
the study of the BEC phenomenon in completely new 
regimes. Polariton condensates differ in several funda- 
mental aspects from the ideal case generally considered 
in textbooks: polariton-polariton interactions are signif- 
icant and the system is far from thermodynamical equi- 
Hbrium. While the separate effect on condensation of 
the interactions and of the non-equilibrium condition is 
already well understood from cither ultracold atom [t'I 
or laser Q theory, not much is yet known about the in- 
terplay of the two effects when simultaneously present. 
In this case, the Bose gas is in fact a quantum degen- 
erate, interacting many-body system whose stationary 
state does not correspond to a thermal equilibrium state, 
but rather originates from a dynamical balance of pump- 
ing and losses [1, [l^, 11 1 

Some striking consequences of the non-equilibrium 
condition on the elementary excitations have been re- 
cently predicted, the propagating sound mode of equi- 
librium condensate being e.g. replaced by a diffusive 
13, Even emarkably, an unexpectedly 
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rich behavior has been observed in recent experiments 
in the spatial and spectral shapes of the condensate de- 
pending on the size of the pump spot [l], Q. While at 
equilibrium, BEC generally occurs at zero momentum 
and the trapping potential is only responsible for the 
/c-space broadening due to finite size [7|, the first ex- 
perimental studies of polariton condensates performed 
with a relatively small pump laser spot showed Bose- 
Einstein condensation into a ring of momentum states 
with a non-zero wave vector il,] ; mutual coherence of dif- 
ferent /c-states was however interferometrically demon- 
strated, which proved that a true condensate was created, 



and not a fragmented one [l^ . Standard condensation 
around A: = was then recovered in later experiments 
using a much wider excitation spot 0,13. So far, most 
of these experimental observations have challenged the- 
oretical understanding [13J: the purpose of the present 
Letter is to propose a complete and unified theoretical 
model able to explain them in a simple and physically 
transparent way. 

Our work is based on a mean-field study of non- 
equilibrium condensates which takes explicitly into ac- 
count the effect of the spatially finite pump spot by 
means of a generalized Gross-Pitaevskii equation (GPE) 
as developed in 11]. This approach can be used to theo- 
retically model non-equilibrium Bose-Einstein condensa- 
tion for any pump laser geometry and for any microcavity 
disorder potential. As long as mean-field is valid, solving 
the GPE provides complete information on the shape of 
the polariton condensate in both real- and momentum- 
space, as well as on its spectral properties under a pulsed 
excitation. The formal similarities with the equations ap- 
pearing in the context of pattern formation in nonlinear 
dynamical systems far from equilibrium 1J|, in particu- 
lar hydrodynamical and nonlinear optical [la| ones 
will be of great utility to obtain a physical understanding 
of the complicate spatial structures that appear as a con- 
sequence of the interplay of inhomogeneity, nonlinearity, 
driving and dissipation. 

At mean-field level, the dynamics of the condensate 
macroscopic wavefunction ip{v) is described by a gener- 
alized Gross-Pitaevskii equation [11] of the form: 



[R[nR{v)] - 7e 



K.t(r) + hg |V;(r)|2 + Vr{v)]4,{v), (1) 



where ~ hujo and m are respectively the minimum 
and the effective mass of the lower polariton branch and 
g > quantifies the strength of repulsive binary interac- 
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tions between condensate polaritons. Whenever needed, 
cavity disorder can be included as an external potential 
term T4xt(r)- Condensate polaritons have a linear loss 
rate 7c, and are continuously replenished by stimulated 
emission from the polariton reservoir created by the non- 
resonant optical pump. At the simplest level, the corre- 
sponding gain rate R[nj^] can be described by a monoton- 
ically growing function of the local density n/j(r) of reser- 
voir polaritons in the so-called bottleneck region ■ At 
the same time, the reservoir produces a mean-field repul- 
sive potential Vfi(r) that can be approximated by the lin- 
ear expression ^^(r) ~ hgRnR{r) + HQ P(r), where P(r) 
is the (spatially dependent) pumping rate and (7_r,5 > 
are phenomenolorical coefficients to be extracted from 
the experiment The GPE equation ([1]) for the con- 
densate has then to be coupled to a rate equation for 



the local condensate wave vector kc(r) is defined as the 
spatial gradient of the condensate phase kc(r) = Vr0(r). 
Stationarity of the solution then imposes that 
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hR(r) = P(r) - 7fl"i?(r) - i?[r^i^(r)] IV^r)]^ 
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polaritons are injected at a rate P{y) and relax at an 
effective rate 'Jr ^ 7c [l^l- Depletion of the reservoir 
density due to the stimulated emission into the conden- 
sate mode is taken into account by the _R[n/j(r)] |'0(r)P 
term. 

In the homogeneous case, i.e. under a uniform pump- 
ing and in the absence of any external potential, the equa- 
tions ([l][2]) admit simple analytical stationary solutions. 
Below the threshold, the condensate density remains zero 
IV'P = 0, while the reservoir one grows linearly with the 
pump intensity, ur = P/^r- At the threshold pump 
intensity P*'', the stimulated emission rate exactly com- 



pensates the losses R[. 
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Provided the spatial variation of the pump profile P(r) 
is smooth enough, one can perform a kind of Local Den- 
sity Approximation (LDA), where the quantum pressure 
term (proportional to V^y^) in ([6]) and the current di- 
vergence term in ([7]) are neglected. Within this LDA ap- 
proach, the condensate density vanishes p{r) = at all 
points r where the pump intensity is below the threshold 
P(r) < Pth- As the oscillation frequency ojc of the con- 
densate is spatially constant, the variation of P across 
the pump spot must be compensated by a corresponding 
spatial variation of the local wave vector kc. 

For a circular pump profile P(r) and in the absence of 
disorder Vext = 0, the condensate frequency iVc is deter- 
mined by the condition that the local condensate wave 
vector vanishes kc = at the center of the spot: Wc— = 
g p{r = 0) -|- gRnR{r = 0) + Q P{r — 0). For standard 
(e.g. gaussian) pump spots with a monotonically decreas- 
ing intensity profile along the radial direction, kc is in 
the outward radial direction and its modulus fcc mono- 
tonically grows in the radial direction, and reaches its 
maximal value at the condensate edge where P(r) = Pth- 
The repulsive interactions create in fact an antitrapping 
potential Vat{r) = Kg p(r) -f hgRURir) + HQ P{r) which 
ballistically accelerates the condensate polaritons away 
from the center and the state equation Q simply de- 
scribes the conservation of the mechanical energy during 
the flow. In order to verify the robustness of these an- 
alytical considerations and assess their validity for the 
typical values of the experimental parameters, extensive 
numerical simulations of the GPE equation ([T]) coupled 
to the reservoir evolution equations ^ have been per- 
formed for a wide range of pump parameters. 

The case of a cw pump with a wide Gaussian spot of 
+ 9\'^gV + 9r''^R + Q P- (3) waist Up = 20 p.m is shown in the four plots in the left col- 
umn of Fig. [1] The maximum height of the antitrapping 
potential Vat is chosen in a way to reproduce the exper- 
imentally observed blue-shift of the emission frequency 
Wc — ^0 of the order of 1 meV 0. As expected, the long- 
time dynamics of the system tends to a steady state with 
a single oscillation frequency Wc and stationary reservoir 
densities. In agreement with the previous analytical dis- 
cussion, the fc-space distribution [Fig. [Ijb)] is contained 
in the k < fcco region (Hk^Q/2m ~ lo^. — uq) delimited by 
the free particle dispersion (vertical dashed lines in the 
figure). Still, the broadening of the condensate momen- 



sate solution with = becomes dynamically unstable. 
Above the threshold, the reservoir density is homoge- 
neous rtfl(r) = Ur and the condensate wavefunction is of 



'00 ■ Inserting this ansatz into 



the motion equation, we find that ur is clamped at the 
threshold value n^, while the condensate density grows 
as l^/ioP = {P — P^'^)hc- The condensate wave vector kc 
remains so far undetermined and stable solutions with ar- 
bitrary values of kc can be found. Finally, the oscillation 
frequency lUc is fixed by the state equation: 



The physics is much richer in the presence of an in- 
homogeneous intensity profile P(r) of the pump. In this 
case, we can still look for stationary solutions of the form: 



V;o(r)e~-^* 
nR{v). 



VW)e 



(5) 



Here, p(r) and (/)(r) are the local density and phase of the 
condensate wavefunction ■0(r) = exp[z0(r)]. The 

condensate frequency lOc is the same at all points, while 
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FIG. 1: (Color online) Numerical results in the absence of 
disorder for a circular excitation pump spot. The polariton 
distribution is shown for a large ap — 20 ^m (a-d) and a 
small (jp = 1 /im (e-h) pump laser spot. Panels (a,e) give 
the (fc, E) emission pattern, (b,f) the polariton distribution 
in reciprocal space, (c,g) in real space and (d,h) the local 
wave vector kc{r). The wave vector kco of a free polariton 
at ujc is indicated by the dotted lines in (a,b,e,f); the dashed 
line in (g) is the analytical approximations to the density tail, 
and the dashed line in (d,h) are the LDA predictions to the 
local wave vector. All quantities in real (momentum) space 
depend only on the radial coordinate r = |r| (fc = |k|). Values 
of parameters used in the simulations: hg — 0.03meV/im^, 
fi-7c = 0.5 meV, hjn = 2meV, hR[nR\ = (O.lmeV/xm^) x na 
hgn^O, g = 0.035^m^ and P/Pth = 2 (left panels), G = 
and P/Pth = 48 (right panels). 



turn distribution due to the ballistic acceleration largely 
exceeds what predicted by mere finite size effects. This in 
stark contrast with the case of equilibrium condensates, 
where the condensate phase is uniform throughout the 
whole cloud and the local wave vector is thus everywhere 
kc = 0|. To confirm our interpretation, we have com- 
pared the numerical result for the local wave vector with 
the LDA prediction. The two are plotted in Fig. HJ^d) 
as respectively a full and a dashed line: the agreement 
between the two is everywhere excellent. 

Although this simple calculation is able to correctly 
account for the main experimental observations, some 
features are still missing. As one can see in Fig. 8. 9 of 
Ref. [iOl , some structure is in fact present on top of the 
broad profile shown in Fig.[T]Ja,b), in particular a narrow 
peak centered at a finite k ^ momentum. To explain 
these features we have performed GPE simulations in- 
cluding a weak disorder potential Vexti'^c) acting on po- 
laritons. Fig[2] gives an example of the outcome for a 
given realization of the disorder potential, constructed 
as an ensemble of circular defects [Figl^l^a)]. The k- 
space emission shown in Figl^fc) is still contained in the 
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FIG. 2: (Color online) Numerical results with disorder, (a) 
Disorder potential (gray scale with ImeV range) , polariton 
distribution in 2D (b) real and (c) momentum space . The 
arrows in (b) indicate the polariton current. The radius of 
the dotted circle in (c) is the free polariton wave vector kco- 



k < kco circle as in the clean system case, but a signif- 
icant speckle-like modulation appears on top of the reg- 
ular background, with a few dominating narrow peaks. 
This phenomenology can be understood in terms of the 
interference of the emission from the different spatial po- 
sitions: in agreement with experimental observations, 
the spatial profile of the condensate [Fig[2lb)] consists 
in fact of several, yet mutually coherent spots connected 
by lower density regions. In particular, the disorder is 
not able to spoil the outward radial flow pattern which, 
despite some local distortions, remains clearly visible in 
the figure and is responsible for the k 7^ position of the 
peaks in the k-space pattern. 

The properties of the polariton condensate are com- 
pletely different if a small excitation spot is used. Even 
though the local density approximation can no longer be 
performed, exact analytical predictions for the conden- 
sate wave function can still be obtained in the region far 
outside the pump spot. In this region, the condensate 
density is very small and polaritons no longer feel the 
repulsive potential, so that the radial part of the GPE 
(H]) reduces to a linear Helmholtz equation of the form 
{h/2m) V^V + {'^c — t^o + «7c/2)V' = 0, whose converging 
solution for r — > 00 has to be chosen. The frequency lOc 
is fixed by the dynamics in the excitation region, and its 
analytical determination goes beyond the scope of the 
present work. The solution to the Helmholtz equation 
has the simple analytic expression il)(r) — H'^^ikcr), in 
terms of the first kind Hankel function iJ^.^' . The com- 
plex wave vector is kc = yj2m{ijjc — ciJo + i'icl'^)/^- From 
the asymptotic expansion of one immediately gets 

to the exponential decay ^-){r — > 00) cx exp(— Kpen'')/\A' 
with a spatial rate Kpen = Im[fcc] ~ 7c/(2wc), and a 
propagation speed Vc = hkc/m, with kc = Re[fcc] — 
^2m{ujc — wo)/fi [llj. These analytical considerations 
are accurately confirmed by a full numerical integration 
of (IT][2]) for a small pump spot CTp = 1 /im, whose re- 
sults are summarized in the plots in the right column of 
Fig. [H As the pump laser spot ap is much smaller than 
the characteristic propagation length k~^^, most of the 
polaritons are ballistically propagating outside the exci- 
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FIG. 3: (fc, E) emission pattern from a simulation with a small 
pump spot (cTp = 1 fj,m) and the following time profile: linear 
rise time of 1.3 ps, plateau and fall time of 2.6 ps. 



tation region. The fc-space emission then concentrates on 
the ring of radius fee! the contribution from the excita- 
tion region provides the weak pedestal at k < kc, while 
almost no emission is present for k > kc- The (fc, E) pat- 
tern is analogously peaked around (±fcco, and shows 
a weaker pedestal on the horizontal segment aX E — hujc 
contained inside the free polariton dispersion. Qualita- 
tive agreement of these results with experimental obser- 
vations in Fig.l of Ref. [l| is excellent: the ring structure 
is recovered as well as the relation between the emission 
frequency lOc and the wave vector k^. 

As a final point of the Letter, it can be interesting to 
investigate the effect of a finite duration of the excitation 
pulse. The result of a numerical simulation of the GPE 
model [i^l is shown in Fig. [3l Although the condensate 
does not have time to reach its stationary state, its basic 
properties remain similar to the ones shown in the right 
column of FigUl The main difference is the appearance in 
the (k, E) emission pattern of additional horizontal lines 
corresponding to new emission frequencies. It is interest- 
ing to note that this feature is in remarkable qualitative 
agreement with the experimental observations shown in 
fig.2b of Ref. [X]. 

In conclusion, we have developed a full theory of the 
spatial and spectral shape of polariton Bose- Einstein con- 
densates which is able to explain all the peculiar features 
observed in recent experiments [l|, |3, US] ■ Significant dif- 
ferences are pointed out with respect to the standard case 
of condensates at equilibrium: because of the repulsive 
interactions, the steady-state is characterized by a signif- 
icant flow in the outward radial direction and polaritons 
can ballistically propagate outside the excitation spot. 
The agreement of the results of our generalized Gross- 
Pitaevskii approach with experiments is promising for 
the application of the present model to more complex 
geometries that may be relevant in the design of polari- 
ton lasers based on suitably patterned structures. From 
the theoretical point of view, the next step will be in- 
clusion of fluctuations around the mean-field, so as to 
obtain quantitative predictions for the condensate coher- 



ence properties and, eventually, for the critical properties 
in the vicinity of the non-equilibrium phase transition. 
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Shortly before submission of this manuscript, the 
preprint [23| has appeared where similar issues are ad- 
dressed from a slightly different perspective. 
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